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ABSTRACT 


The equations of motion of an arbitrary flexible body in orbit 
are derived. The model includes the effects of gravity with all its 
higher harmonics. As a specific example, the motion of a long, slender, 
uniform beam in circular orbit is modelled. The example cons5.ders 
both the inplane and three dimensional motion of the beam in orbit. 

In the case of planar morion with only flexural vibrations, the pitch 
motion is not influenced by the elastic motion of the beam. For large 
values of the square of the ratio of the structural modal frequency to 
the orbital angular rate the elastic motion is decoupled from the pitch 
motion. However, for snail values of this ratio and small amplitude 
pitch motion, the elastic motion is giverned by a Hill's 3-term equation. 
Numerical simulation of this equation indicates the possibilities of 
instability for very low values of the square of the rario of the modal 
frequency to the orbit angular rate. Also numerical simulatiors of the 
first order non-linear equations of motion for a long flexible beam 
in orbit have been performed. The effect of varying the initial con- 
ditions and the number of modes has been demonstrated. 
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1. IbTRODUCTION 


This report presents the development of the equations of motion 

of an arbitrary flexible body in orbit and its specialisation to a 

12 3 

long, slender, uniform beam in circular orbit. In the literature, ’ ’ 
many models of a beam satellite in orbit are presented. The model 
considered by Pringle consists of a massive central rigid body to 
which a massless elastic beam is rigidly attached. The other end of 
the beam carries a tip mass. For beam satellites with lengths of the 
order of 100 meters or more and mass distributed throughout the length 
this cannot be used. Ashley’s model ’ of a beam satellite is based on 
a continuum approach. He arrives at a set of partial differential equa- 
tions for the beam model through energy methods. However, partial dif- 
ferential equations are not as convenient for numerical simulation. 

Ashley arrives at the conclusion that it is impossible to excite flexural 
vibrations of a beam directly through gravity gradient. The same con- 
clusion is arrived at in section 5.1.1 of this report. 

The development of the beam model presented in this report follows 
as a specialization of the equations of motion of an arbitrary flexible 
body in orbit. Also, the model presented here can be quite easily adopted 
for numerical simulation and also for inclusion of various control laws. 

The development of equations of motion for an arbitrary flexible 

4 

body presented in this report essentially follows that of Santini. 


However, the distinguishing features of the present development frcm 
that of Santini are (a) extensive use of the techniques of vector 
calculus (b) introduction of the orbit fixed reference frame as an 
intermediate frame between the lcoal intrinsic frame at the origin 
of the body axes system and the body axes frame. Hence, the Euler 
angle rotations used in the present report are different from those of 
Santini. 4 

The equations of motion presented in this report are obtained by 
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the Galerkin integration ’ of the equations ox motion of a generic 
point in the body. Hie motion of the generic point is assumed to be 
described by the superposition of rigid body motion plus a combination 
of the structural modes. 

Section 2 of this report describes the various co-ordinate reference 
frames employed in the development of the equations, and the transformation 
relations between these reference frames. The detailed derivations of 
the transformation relations are presented in Appendix - I. 

Section 3 presents an expression for the gravitational force on 
a generic element in the body. The expression also includes the higher 
harmonics of the earth's gravitational field. The details of this 
development are presented in Appendix - II. 

In Section 4, the development of the equations of motion of an 
arbitrary flexible body in orbit is presented. The expansions of the 
vector expressions in equations (28), (31), (32), (43), (44), (4S), and 
(4S) are presented in Appendix - III. 


Specialisations of the equations of notion in section 4 to a long, 
slender, uniform bean are presented in section 5. Sections 5.1 and 
5.2 discuss the planar, and the three dimensional notion of the bean 
in orbit, respectively. In Appendix - IV, the expressions for the 
natural mode shapes, frequencies and the nodal mass for a free-free 
uniform beam are presented. 

In Section 6, the numerical solutions to the response of the planar 

motion of a beam in circular orbit, undergoing only flexural vibrations, 

are presented. The responses indicate the possibilities of instability 

2 

at very low values of (w n /o) c ) . 

The final section of the report, Section 7, outlines the conclusions 
based on the numerical results in Section 6. 



2 . CO-ORDINATE FRAMES 




The following co-ordinate frames are used in the development of 
the equations of motion. (See Fig. 1); 


: O'XYZ Inertial reference frame. 

0 f Z along the earth's spin axis. 

O'X along the line of the ascending node. 
O'Y perpendicular to O'X and O'Z. 



Pi^ijig bocal intrinsic frame at a generic point, P. 

Pi^ along the radius vector from O' to P. 

Pi 2 perpendicular to Pi^ in the plane of ZO'P. 
Pi^ perpendicular to Pi^ and Pij* 

Ul l~2“3 local intrinsic frame at 0. 


?2 : OXqYqZq Orbit fixed reference frame. 

OXp along the local vertical 

OYp along the orbit normal and in the negative direction 
of u the orbit angular momentum vector. 

OZp perpendicular to 0X^ and OY^. 

t 3 ; OXYZ Principal axes of the body. 


The above reference frames are related to each other as follows: 
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0 c.x 

S X 



_2 ~ S X 

c JC 



C<J)C0 

CS(f)C^ + O0S0S1^) 

(stJjsiJ* - C0S0C^) 

-stfice 

Cc0c^ - s<J>s0sij;) 

(c4>s^ + s<}>s0a(ij 

s6 


-C0stj/ 

e0c^ 


( 2 ) 


(3) 


and s,c represent sine and cosine functions respectively. The detailed 

development of these transforrrations are presented in Appendix - I. 

The body angular velocity components (w , w . , w ) and the Euler 

x y z 

angular rates ($>, 0, $ ) are related as follows: 


to = 0s<{> + ^ccj»c0 - 0 ) (S<J>C^ + c<|>s0s^) 

X' c 

Uy - 0C<J> - ^S<|)C0 - (1> C Cc^c^ - S0B@3lp) 

w = ^s© + <£ + w i c0si/j (4) 

w C 



3. GRAVITATION 


In this section an expression for the gravitational force per 
unit mass, expressed in the tody fixed frame (t^) is presented. 

In deriving this expression, it is assumed that |r|/p« 1, where r 
is the position vector of an arbitrary point in the body with respect 
to 0 and p is the distance of the same point from 0 ’ , 

We can write the gravitational* 4 potential in the most general 
form as ; 


where, 


v (p, n, = 


\>a* 


+ va I 

S-l 


«. (i-J 


S+l 


n (n,w) (5) 

S 


K. = cos A 
s so T so 


fl (n,w) = l [F c Cm) (n) cos + <J> e J]/K c 
m=0 s 


'sm 


P Cn) is the associated Legendre function of 
order s, 

K go and # are constants to be given experimentally through 
geodetic satellite techniques. 

The gravitational force per unit mass at the origin of the body 
axes, 0, in the t^( 0) franeis F^ - V V|q. 

For a point at a distance r from 0, neglecting snail quantities of 
the order of the gravity force in the x q frame is given by 

P 0 " “ 


F « Fq + M r 


C10> 



where 


f = Gravitational force at 0 expressed in the principal body 


axes frame (tg) 


<n) UaT 

M = p 3 


[ M (0) + I 

K s (— ) S 

M (s) ] 

s = 

i s V p / 






2 2 

3c (J>c 6-1 

2 

— 3s<Jjc<J>c 0 

3c<f>c6s0 

2 

-3s4>g<{>c"0 

5 2 

3s $c e-i 

-3s<t>c0s0 

3c<{>c0s0 

-3s^c0s0 

3s 2 B-1 

— 


' 

.. 2 

(c) m 

(s) T 
To T« B vs; T 9 

_ T 
To 


( 6 ) 


The matrix B^ s ' is given by eqn. (11-12) in Appendix - II. It can 


be observed that and are synrnetric matrices. 


The expression for the gravitational force in the form given in eqn. (6! 
is used in the development of the equations of motion presented in the 
next section.. 


4. EQUATIONS OF MOTION 


In this section the equations of rotational motion and elastic 
motion of an arbitrary flexible body in orbit are presented. The 
body is assumed to be subjected to small amplitude elastic displace- 
ments, q, which are transformed to elastic forces by the linear 
operator L. 

Consider an elemental mass, dm, whose instantaneous position 
vector from the center of mass of the body is r (Fig. 2). 

The equation of motion of dm can be written as: 

a dm * L (q) + F dm + F (7) 

where 

a = Inertial acceleration of dm 
L(q) = Elastic forces acting on dm 

¥ = Gravitational force per unit mass 
E " External forces acting on dm 
q * Elastic displacement 

The above vector equaton can be written in the body fixed reference 

frame (tg) as, 

.. • • 

[a^ + ?+ 2 Jxr + Jxr + ux (wx r)] dm 

= L (q) + F dm + E (8) 

It is important to note that r and r are the velocity and acceleration 
of dm respectively as seen from the body fixed reference frame, tg. All 
the vectors in the above equation must be expressed in the body frame 

V 
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We can write the instantaneous position vector r of dm as 


r = r n + q 


( 9 ) 


where 

Tq is the position vector of dm with respect to 0 in the 
undeformed state, 1 q is the elastic displacement of dm. 

Hence 

• • *» »« 

r = q and F = q (10) 


For small amplitude elastic displacements, we can write the elastic 
displacement, q, as a superposition of the various modal contributions 
according to 

CO 

q = Z A(t) ^ n) (? n ) (11) 

n=l n u 


where 


^ n) C V = 4 ’x <n) 


i + 4 


Cn) 


j + 


(n) 


T^n) 

<p 


A(t) = Modal amplitude 

(r^) is Idle mode shape associated with the natural frequency 


w and satisfies the following orthogonality condition, 


Also 


i i (m) 
vol 


— (n) 

rda = <5 M 

mn n 


L C* 00 ) 


2 

-u $ dm 


( 12 ) 


(13) 


Further, if the body is unconstrained, the elastic modes must be 
orthogonal to the rigid body modes ie, 


f dm = 0 (14) 

vol 


- 9 - 



( 15 ) 


/ r n x $^ n; dm = 0 
vol u 


If the body is constrained against translation and rotation at the 

undeformed center of mass, the corresponding inodes are called "fixed 

modes." For fixed modes the orthogonality conditions (14) and (15) do 

not hold. It should be noted that, for the case of fixed modes, the origin 

of the body frame, 0, no longer coincides with the center of mass in 

the deformed state. Hence, only for free modes, / r dm = 0. However 

vol 

for fixed modes S r dm = / q dm ^ 0. 
vol vol 

4.1 Equations of rotational motion 

The equations of rotational motion of the body are obtained by the 
following operation 


/ r x Eqn. (8) 
vol 

i*®* .. . « 

/ r x [a_ + r + 2oxr + wxr + wx (w x r)] dm 

vol OT 

- / r x [ L(q) + F + e] dm (16) 

vol dm 

where e is the external force per unit mass. 

The various terms in Eqn (16) can now be evaluated using the techniques 

of vector calculus. Assuming [q| « 1, only the first order terms in 

ir| 

q are retained in the following expansions. 


/ 

vol 


r x a dm = 

cm 


; 

vol 


q dm x 
* cm 


(17) 
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( 18 ) 


*• i» »» 

/ r x ? dm = S (F.+q) x o dm < / r Q x q dm 

vol v u v u 


/ r x 2 (u x r) dm * 2 /(r n +q) x (w x q) dm 
vol V u 

ft 

= 2 / r n x (w x q) dm 
v u 

% ft 

/rx (wxr) dm = / (x\.+q) x [« x <r 0 +q))dni 
vol V 

ft ft 

= / r Q x (to x r 0 )dm + / r Q x (w x q)dm 
v u u v u 


♦ / q x (to x r Q )dm 
v u 

/ r x (to x (to x r))dm = / [(r * to x r)to - (r » w) 

vol v 

= - / (r • to) (to x r)dm 
v 

= - / {(r^+qV to) (to x (rQ+q))dm 


/(r Q * to) (w x r Q )dm 


- / t (r 0 • S) (to x q) + (q • to) (to x 
v 1 


/ r x L (q)dn\ = / r x L (q)dm 
vol V 0 

W /j. V 

= - I to n " ^ x ^ ; dm (using eqn. 
n=l 


(19) 


( 20 ) 
(to x r)]dm 


r^dm (21) 


(13)) 


0 (for free modes * See Eqn. (15)), 


( 22 ) 


/ F x F dm » / r x (f n + M r)dm 
vol v u 


= / q dm x f n + / (r Q +q) x M (r Q +q)dra 
v u v u u 


/ q dm x f 0 + / ?q x M r Q dm 
v v 


+ / [r Q x M q + q x M rJdm 
v 


(23) 


/ r x e dm = C 
vol 

It can be easily shown that, 

* » A t a ^ A 

/ r Q x (u x r Q ) dm = w x i + J y u) y j + J 2 w z k (24) 

A A 

and / (? 0 • w) (w x r ft ) dm = ui xr ui 9 i + (J^ - J v ) w, j 


y ~z' y ~z 


X Z X 


+ <J x - V “x “y k 


(25) 


where, J , J , J are the principal moments of inertia of the body in 
x y z 

the undeforroed state. 

After substitution of the values for the integrals in Eqn. (16) and 
rearrangement of terms, one can obtain the following form for the equations 
of rotational motion. 


DO 00 CD 

R + 2 + C + I = G_ + 2 

n=l n=l * n=l 


(26) 


where. 


5 * {J x ^ . <J z - J y ) Wy Mj ,} i 
+{J y «y * (J x ' V “z V 1 

* fJ Z i + (J y - J x > U X “y } k 


( 27 ) 


I tf n) = / [r n x q + 2r n x(uxq) 

VGl U u 


+ r Q x (w x q) 


+ q x (w x r Q ) - (? c • w) (US’ x q) 

- (q • w) Cw x F 0 )J dm (28) 


G a / r x e dm 
vol 


(29) 


I 

n=l 


D <n> = / 


q dm x (a - F ) + Z 
vol n=l 


CO 

n 


f F 0 x ^"'d m (30) 


Sr = f Fq x M Fg dm (31) 

vol u 

Z ^ n) = / [r n x M q + q x M ? n ] dm (32) 

n=l vol u u 


The expansions of the above vector expressions are presented in 
Appendix - III. 

The significance of the various terms in Eqn. (26) will now be briefly 
discussed. The terms, ($<»> , reflect the inertia torque associated with 
the elastic deformations. The term, G^, corresponds to the gravitational 
torque on a rigid body. The terms, , correspond to the gravitational 
torque due to the elastic deformations. The terms, Ty , account for 
the difference in position between the actual center of mass and the 
center of mass of the undeformed body. For the case of free modes, 

D* n) = 0. 
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4.2. Generic mode equations 


The generic node equation is obtained by the following 
operation: 

/ ^ n) *C£qru (8)3 

vol 


le 


/ F • [a + r * 2 w x r + w x r + w x (to x r\>] dm 
vol ^ 


= / ?< n) . [tt|L + r + si 


dm 


(33) 


vol 


With the use of equations (6), (9), (10), (11), (12) and (13), 
the various terms appearing in Eqn. (33) can now be expanded as follows: 


/ a dm = / T dm * a 

vol 0111 V m 

S ♦ r dm • f * q dm 
vol v 


= l \S ? (n) - ^ m) 


m 


dm 


(34) 


/ $ n »(2toxr)dm :: 2/ w x q) dm 


vol 


(35) 

(36) 


a q £ 

/ (4*^‘ to x r) dm = / C^ n ^* ST x r 0 ) dm + / (<F^ n to x q) dm (37) 
vol v v 


/ w x (to x r) dm = S ♦ to x (to x r n ) dm 

vol v u 


+ / $ n ^« u x (w x q) dm 


v 


(38) 


V 


/ 5<n> 


1 vol 


• L Cq> = - r U 2 A / ^ n) . *<»> dm 
m=l m v 


= - u) A M 
n n n 


/ ^ n) * f dm = / ^ n) dm • f_ + / $ (n) . m r 
vol v u v 


dm 


xCn) 


= 1 $V1W dm • f _ + / ? (n) . M r A dm 


v 


‘0 


v 


+ / m q dm 


/ ^ n) * e dm = E 
vol n 


(39) 


(40) 

(41) 


A^ter substitution of the values for the various integrals in Eqn. (33) 

and rearrangement of the terms, the generic mode equation is obtained in 
the following form. 


a*» 2 a *% 


+ — 1 < s , mn =F^-[G n + E n ♦ d;: 


where 


<9 n = ^ f^ n) * u» x r Q + ^ n) . W x (wx r Q )3 dm 


(42) 


(43) 


\*hn a ! [2* <n) - u X ? + 5r< n5 
m=l v 

• 

. w x q + T (n) . u x (w x q)] dm 

g n = f ^ n) - M r Q dm 

(45) 

00 

£ g= / $ Cn) • M q dm 
m=l ^ nn v 

(46) 

E 5 / e dm 

V 

(47) 

D' = / * <n) dm • (a - 7 n ) 
n cm o 

(48) 

M n = / T*™* dm 

(49) 


(44) 
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The expanded ferns of the above expressions are presented in 
Appendix III. 

The significance cf the various terms in Eqn. (42) will now be 
briefly discussed. The terro^ corresponds to the forcing tern due 
to rigid body motion. is the forcing term due to the elastic 
motion in the m mode. The term, ff , represents the gravitational 
force acting on the mode due to the rigid body motion. The terms 
g correspond to the gravitational force acting on the n ‘ mode, due 
to the elastic motion in the m tn mode. E n is the component of the ex- 
ternal force acting on the n^ 1 mode. IT is the term corresponding to 
the displacement of the center of mass from the point 0. 

In the next section, an application of the equations of motion 
developed in this section is presented. 
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5. APPLICATIONS 


In this section , we consider the application of the equations 
of motion presented in the previous section, to the specific example 
of a beam in circular orbit shown in Fig. 3. 

In section, 5.1 it is assumed that the motion of the beam is 
restricted to the orbit plane . In section 5.2, the general case 
of the three dimensional beam is considered. 


5.1. Planar motion of a long, slender, uniform beam in circular orbit 

Since the beam motion is restricted to the orbit plane, the 
yaw and the roll angles vanish, i.e. 

tp (t) = $ (t) - 0 (50) 


Also, the elastic deformation out of the orbit plane is assumed 
(r ) 

to be zero, ie, 4> y - 0. 

For unconstrained structures. 


d n) = OS p; > o and - h£> 


(51) 


= / £ 4 n) dm (a, S = x,y,t) 

aS vol a S 


and (£,£,£) are coordinates of dm in the unde formed state measured 
x y z 

in the body frame. 

In the absence of out-of -plane deformations, i.e. 4^,^ = 0* we 
can deduce that, 

(n) . H (n) _ w (n) . H (n) _ „(n) _ n 
V ' K XV - V ‘ H 2V - H V2 - 0 


( 52 ) 


Noting that, for a bean (see Pig. 3) with uniform cross section. 


* * r Jr j ? 


we can further deduce 

..(n) . , 

n - ' 

- ,Cn) 

cx 


vol 


p/ n ' = I 

,.r .<n) 

VU._ <f w 

z ~ vol 


where, y is 

rile densi* 


^x 5 C S< 


(5>d£ 

X x 




= 0 = t 


_ -.Cn) 


xc 


;S_ = 0 


( 55 ) 


Since, bv definition. 


L (nr ‘ ? = r g fe) ti Vi) dm (a, 5 = x, v. -) 


“c\S 


r 

vol 


. . , - _ , , . , , _ . . (a) . 

tn tne absence or oub-o.r-p.tine exastac d.ercrraatxon, x.e.. * - o. -t. c^n 

be easily shown that 


. ten) _ , (inn) _ . (m) _ . (nr.) _ . (mn) 
bvV S*x ‘ ~Vc ' tev " Hv 


(55) 
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tne results rresentec 


c_ tie above results and 
Appendix - III, the pitch and tne generic node equations for the pre- 
sent case are obtained as: 
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The generic mode equatiore: 
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+ a 2 A + %. + 1 

r n u 
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(58) 


where, 

% - - (5 - “/ 

■ I 2 \ - “c> + \ «] «£ 1> - g°> 


- A (8 - w ) 2 (L (mn) + L (ran) ) 

m C XX 22 

^ = “n I *»! > 

V = \ % ^ + «U CL “ + + M 33 L “l 


We will now consider several specific cases of the long slender beam 
assumed to be subjected to rotations and deformations only within the 
orbital plane. 

5.1.1. The case of no longitudinal vibrations: i.e., <j/ n ^ = 0. 


For this case it can be seen that, 


H. 


(n) _ L (mn) _ L (mn) _ L (mn) _ Q 


xx 


xz 


zx 


XX 


(59) 


(n) 


By choosing t}> to represent the eigen-modes of bending vibrations 
z 


of a free- free beam, we obtain, 

= fi M 
22 mn n 


(60) 


where, 6 is the krone cker delta and M is the generalized mass of the 
mn n 

beam in its mode. An expression for the value of for free-free 
beams is presented in Appendix - IV. 


Substituting the above results in equations (57) and (58), one 
obtains the pitch equation and generic mode equation in the following 


form. 

8 + M,. + 1 * 0 
31 J y 

Generic mode equation: 

*n * K 2 - m 33 - ■ r 

n 


(61) 


(62) 


where, 


M„ = ^| Cs0 c0 + Z 
31 p 3 


K 


s~l 


( a.) s M (s) l 

V n n i 


«33 = ' 3 s2 6 * 1 + j 2 K s ( | )S «H ) 


Assumptions of a spherically symmetric gravitational field and 
circular orbit result in the following simplification of and M^, 



M„, = t w 2 sin 20 
31 l c 

M- 0 = to 2 (sin 2 9-1) 

3o c 

Hence, equations (61) and (62) simplify to, 


6 + I- w 2 sin 20 + C v =0 (63) 

2 c 

A n + (o) n 2 - o) c 2 (3 sin 2 6 - 1) - (6 - w c ) 2 ] A^ * ^ (64) 

Equation (63) is just the rigid body pitching equation. Thus from 
Eqn. (63), one can conclude that there is no influence of the elastic motion 
on the rigid body motion, under the assumptions of the present case. The 
same conclusions were arrived at by Ashley for the case of a thin beam 
undergoing flexural deformations in the orbit plane. 
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If the rigid body pitch oscillations are small ie, 6 « 1, 
equations (63) and (64) further simplify to (assuming no external 
disturbances) » 

it ^ 

8 + 3 u Z 6 = O (65) 

c 

A^ + [ % 2 - ©2 + 2 w c 6 ] « 0 ( 66 ) 

Without loss of generality, one can assume the solution of equation 


(65) as, 

0 = c sin (/3 us c t + y ) (67) 

where, c is the amplitude of the pitch motion. 


Substitution of Eqn. (67) into Eqn. (66) yields. 

An + [w 2 - 3 w 2 c 2 cos 2 (/3 wt + Y ) 

Ti k n c c 


+ 2/3 w c 2 c cos(/3 ui c t + 1 )] = 0 


( 68 ) 


With the introduction of dinars ionless variables , T “ — U'o* us c t + Y ) 
and z n = A /% where £ is the length of the beam, one can write. 






(69) 


After substituting Eqn. (69) into Eqn. (68), one obtains the 
generic mode equation in the following non-dimensional form, 

2 


d Oh 


+ i [(% - 3 c 2 cos‘2t + 2 /3" c cos 2t] z = 0 

- 3 1 n 


dr 


(70) 


7 1 

Using the trigonometric identity, cos’ I t = ^ (cos 4 t + 1), Eqn. (70) 
can be re-written as, 


7 

d z n 

— 

dT" 



3c 2 ) 


!»Tc oos 2t - ^ c’ cos 4t] z - 0 (71) 

2 n 


av 


Equation (71) is in the form of "Hill’s 3 - term equation" or 

7 

"Whittaker's equation"' For small pitch amplitudes ie. , c « 1, the 

above equation can be further approximated by the Mathieu equation 

z. 

dz n 

— =- + (5 + e cos 2x) z n = 0 (72) 

dt 2 n 


where, 


6 

3 0 ) 




e 


81^3 c 
3 


= 4.62 c 


Figure 3 shows the Mathieu stability diagram 


8 


with e and 6 as 


parameters , It can be seen from the Mathieu chart that for the values 

ajj-,2 , 

of 6 around unitv i.e., — - 5 - = tf ( 1 ) the system may enter the region of 

«cT 2 

instability. However for large values of (i%/w c ) the co-efficient of 
2 n inEqn. (72) is dominated by . Hence, in the high frequency 

range, the elastic modes are essentially governed by the following simple 
equation : 


A 

dT 2 


+ 5- e2Sk >■ 

3 Wc 


= 0 


(^ » 1 ) 

Wc 


2 

Thus, it can be concluded that for beams with (u^Aoq) » 1 
(n = 1 , 2 , ... »), the elastic motion and the rigid body pitching motion 
are completely decoupled from each other. Kcwever, if (u^/o^) = 9 ( 1 ) , 

(highly flexible beams) one has to consider Eqn,(71) to study the elastic 


motion, and thus the elastic motion is coupled with the rigid body pitching 


motion. 


In Section 6, numerical solutions of Eqn. (71) are presented, for 
some typical values of w n /w c . 

5.1.2. The case of no flexural vibrations: (d> ' =0) 

z 


In this Section, we specialize equations (57) and (SB) by restricting 
the beam to have no bending or flexural vibrations, but allow for the 
possibility of longitudi n al vibrations. With the above assumption and 
the assumption of circular orbit in a spherically symmetric gravitational 
field, the following simplifications result. 



= L 


(mn) _ .(mm) 
zx “ xz 


= 0 


M-, = u 2 s0 00 

31 C 

Mn = ** C3 ° 2e " 15 


(73) 

(74) 

(75) 


Substitution of the above results in equations (57) and (58) leads 


to, (assuming C y = 0), 


J y ' + 2 Jl - W C 5 + \ e '] ** * t J y + n Ii H £ )j ^ 2 sln “ = 0 


(n) 

A + w 2 A - [(0 - u ) 2 + a) 2 (3 c 2 0-l)](-~£ + A ) = r r 
n n n c c J M n n M 

n 


(76) 

(77) 


Equations (76) and (77) are consistent with Ashley's equations (A- 8) 
and (A-9 ) in Ref . 2 , which were derived from an energy approach . 

Ashley has further shown in Ref. 1, that if the beam is spinning at 
a high enough spin rate, ie. w c /ft « 1, where ft here represents the 
spin angular velocity, the response of Eqn. (77) for only the first node, 
consists of steady stretching on which can be superimposed an oscillation 
with the physical frequency of (w^ + 2.93ft 2 ) 15 where, is the funda- 
mental frequency of the beam in the longitudinal vibrations assuming 
no spin. 
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Ashley has also investigated another special case of equation (77) 

by assuming 0(t) = constant. Physically this implies that the beam 

is forced to naintain a constant orientation with respect to the 

local vertical. In this case, it can be easily seen that Eqn. (77) 

reduces to that of an harmonic oscillator with the physical frequency 
2 2 2 2 

of K - u) c + (1 - 3c 0)] . The third term in this frequency 

expression is due to the gravity- gradient effect. Thus, the gravity - 
gradient contribution to the frequency changes sign when c § = |. 

For beams with = ©^(u^), the possibility of buckling instability 
at 6 = 0 is also evident from the frequency expression. 

5.2. Three dimensional motion of a long, slender, uniform beam in 
circular orbit: (assuming yaw to be zero, ie. = 0). 


This section presents the equations of three dimensional motion 
of a beam in circular orbit. Seme of the following results which were 
obtained in Section 5.1 also hold for the present case, ie. 


= Os D' n = 0 and = H<£> 


(51) 


Tj(n) u (n) „(n) »,(n) »r(n) »j(n) _ rj(n) »j(n) _ « 

H - ’ H zx - H xy - V - V - zz - yz * H zy * 0 


xz 


Also, J = J and J - J = J for long, slender beams, 
y z y x x 


Hence, the pitch equation is, 

* ,(n) 


J (d - J oj uj + E 0 X 1 1 * + C — S-, 4 E G 

y y y z x n-l ^ y % n-1 y 


(n) 


(78) 


(79) 


- 24 - 








where, 


ai =(6 - to ) s4> i oi- = (6 - to) g<{> ; to = } 
x c y c z 


t^ n> = 2 [A n ui y + A n ^ - 2 a z o> x )]H. 
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xx 


% = - J z M 31 

G (n) = - 2M,. A H^’ ) 
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The roll equation is, 


where. 


J w + J w oi + Z Q (n) + C = G* + E G 
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The generic mode equation is, 

A + to 2 A +5o-..+ 1 
ti n n 1^ 

where, 

d?_ = - (oi 2 + to 2) H fn ^ 

mi ' v z - xx 


— i - — — * 
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5.2.1 The case of no longitudinal vibrations: (<tj n ^ = 0). 

In this case, the following further simplifications result, i.e. 
H Cn) = L (mn) . L (nn> , ,0m) . „ (82 ) 

XX 2X XZ XX 

Furthermore, if we choose the eigen modes of bending vibrations of 
a free- free beam to be * , then, 

f «*w% ^ 

(83) 


L <mn) = 8 M . 
zz mn n 


Assumption of a spherically symmetric gravity field results in. 


= - 3 u) Q 2 s$ c<t> c 2 0; M 22 = (3 s 2 <J> c 2 0 - l)^ 2 
= - 3 w c 2 s4> c© s0 ; ‘ C3 s 2 0 - l)^ 2 


(84) 


M 31 = 3 V 


c<t> C0 S0 

With the above results substituted in equations (79), (80) and (81), 
one obtains the following equations of : 

Pitch: 

(85) 

( 86 ) 
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Generic mode: 
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where, 


oj = (6 - to ) s4> 

*V V 

= (9 - to.,) c<J> 
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The assumption of snail amplitude pitch and roll, i.e. 8 « 1 and 
4> « 1, yields the following set of linearized equations for pitch and 
roll, after neglecting the terms containing the products and powers of 
<J>, }, 8 aid 8. 


Pitch: 


6 + 3 a) * 0 + 
c 


£ 


( 88 ) 


Roll: 


2 

<f> + 3w * + •=— = 0 
c 


The generic mode equation is. 


En 


A + [oj 4 + 2 <f> (8 - to) - (9 - o ) 4 - <J> 4 + 2 c> 4 ]A = M 
n 1 n c c c J n M. 


n 


2 — 

[ 

n 


(89) 


(90) 


5.2.2 The case of no flexural vibrations (<}/ n ^ = = 0): 

y > 

In this case, the assumption of no flexural vibration results in. 

Cm) = L (mn) a jCnn) s 0 (91) 

zz zx xz 

Also, the assunption of to be the eigen-modes of longitudinal 
vibration of the beam leads to the result, 

^ ••• "> 02 ) 


mn n 


As in Section 5.2.1, the assumption of a spherically symmetric gravity 
field, yeilds 

= (3 c 2 <J> c 2 0 - l)w c 4 ; M^ 2 = - 3 s<> c<J> c 2 0 

M 22 = (3 s 2 * c 2 $ - l)u> c 2 : M 23 = - 3 u> c 2 s<$> c9 s6 

M 23 S 3 c6 ^ » M 33 = C3 &2q " 1)a) c 2 


(93) 


For this case , the pitch and roll equations rei/ain the same as 
equations (79) and (80). The generic mode Eqn. , (81), simplifies to the 
following form: 


A + a) 2 A + %- 
n n n hi 
n 


+ t n=l 2 4 <S " 


+ J, E “ + E n 5 

m=l 


(94) 


where, 

% = - s 2 + “ 2 2> C 

= “ A™ [u> } + U 2 ] M 6 
“mn m y z n mn 


s mn = A ta M U M n d W 1 . 


6. NUMERICAL RESULTS 


In this section, the numerical solutions of equation (71), fora 

2 

few typical values of (a> n /w c ) shewn in Table -1, are presented. Cases 
1-5 correspond to very flexible beams. Eqn. (71) was integrated, 
using a Runge-Kutta fourth order method with variable step size, on a 
Nova 840 digital computer. As an initial guess a step size of approxi- 
mately VsQOth of the orbital period was chosen for all the computer runs. 
A pitch amplitude of 0.2 radians was assumed, which is the upper limit 
for 0 for which the approximation sin 6 s' 6 is still valid. 

2 

The responses shown in Figs. 5 and 6 for the value of (o^/u^) =1.0 
indicate the instability of the system at very low natural frequencies. 
Point 1 corresponding to this case on the Mathieu stability diagram 
in Fig. 4 also indicates this instability. The ether points on the 
Mathieu stability diagram correspond to other cases shown in Table 1 
and are referred to by the appropriate case number. For the other 
values of Co> n /u> c ) considered, the beam response is sinusoidal to a 
very close approximation » see Figs. 7, 8, 9 and 10. Phase plane plots 
of cases 1, 2, 3, 4 are shown in Figs. 11, 12, 13, and Irrespectively. 

It is to be noted however, that when the pitch amplitude exceeds 
0.2 radians, one has to use equations (63) and (64) for simulation. 

The development of a general computer program which treats the combined 
pitch and generic modal equations will be discussed in Section 6.1. 
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Case 

No. 




Assigned data 


Conrnents 

(to /to 
K n c 

c* 1 * 

z n <0) 

dz n (0)/dt 

*'y 

1.0 

0.2 

0.5 x 10“ 4 

0.0 

Range of t: 0 to 11 

1.0 

0.05 

0.5 x Ml" 4 

0.0 

Range of t: 0 to 11 

2.0 

0.2 

0.5 x I0 -4 

0.0 

Range of t: 0 to 11 

5.0 

0.2 

0.5 x 10 4 

0.0 

Range of x: 0 to 5.5 

10.0 

0.2 

0.5 x 10" 4 

0.0 

Range of x: 0 to 5.5 

3200 

0.2 

0.5 x 10“ 4 

0.0 

Range of x: 0 to 1.1 


* This case corresponds to a beam with a fundamental natural frequency of ops and 

moving in a circular orbit of 250 n. miles altitude. 

* c - Pitch amplitude in radians . 


Table 1: Data assumed for numerical simulation 



6,1 The Flexbeam Computer Program 


A. Development of Computer Program 

The first order non-linear equations of motion for the long- 
flexible beam in orbit as developed in section 5, Eqs. C57) and (58) , 
were coded for numerical simulation using the Nova 840 digital computer 


system. The program, Flexbeam, consists of nine subprograms whose 


names and functions are: 


FBSET 


FBHST 


FBLST 

FBDIF 


RKSCL 


(a) set up initial state, (b) an input quantity. Precision, 

-5 

dependent upon the desired precision (range 1,0 to 10 " ) ; the 
size of this input parameter is inversely proportional to the 
desired accuracy (and also the number of iterations required 
per computational time step) . 
set where n = 1, 2, . M, * 10 modes 

set where r = 1, 2, K, and n = 1, 2, . M. 

^ j 

set differential equations which include primary and secondary 
functions 

(a) set full scale values of the initial conditions: 4,, <?p, 

* 

An^, and An^, (b) set the number of ordinary differential 
equations to be evaluated, (c) set the bounds on the maximum 
number of iterations (11) that the Runge-Kuita subroutine is 
allowed per time step in oraer to fulfill the desired precision 
specified in RESET. If more than this number of iterations 
would be required the Runge-Kutta subroutine is automatically 
terminated. 


RKGS 


In this case, the precision parameter or full-scale values 
would have to be adjusted. 

Uses the Runge-Kurta method to obtain an approximate solution 
of the system of first order differential equations, given 
initial conditions. 


FSPIF - 


FBSIF - 


FBOUT - 


calculation of primary' intermediate functions, i.e. 
Mu, ^ 13 ’ ^2* cos $> etc. 


calculation of secondary intermediate functions, i.e. 

V Sr’ *m Sm 

sets the scale nagnitudes of the output variables ij>, <}>, A^, and 
A^, equal to the maximum expected amplitudes for use in the plot 
routine. 


The flew diagram, Fig. 15, depicts the various subprograms aid how 
input information, as well as the resulting output of each subprogram, 
is passed through the Flexbeam main program. 


B. Flexbeam Usage 

Flexbeam is a computer simulation of the dynamics of a flexible 
beam in a circular orbit about the earth. The programming language 
used in the min program is Fortran Five. This language allows the 
results to be obtained faster than the same results obtained, using 
Fortran Four. (The Flexbeam out-put program is coded in Fortran Four.) 
Flexbeam, contains three principal components : Flexbeam Main Program, 

(refer to Fig. 15), Flexbeam Input Bata and Flexbeam Output Program. 
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Inside the main program the subroutine FBINP is associated 
with Flexbeam Input Data. The purpose of this subroutine is to 
specify the order and format with which the input data is read into 
the computer. Block data, also inside the nein program, sets the 
constant parameters, i.e. gravitational constant, NU, orbit frequency, 

WD, external forces, QCP, C2, and terms involving the higher harmonics 
of the earth’s gravitational field, J 2 and J3, to their full-scale 
values; also constant parameters associated with each mode are set, 
i.e. the modal frequency, DW. The number of modes CM) and the number 
of differential equations being considered, NORD, are also specified here. 

Flexbeam Input Data is the component program that allows the main 
program to be executed. The first input of the first data card denotes 
the total time interval over which the numerical integration is to 
be performed, e.g. 0.0 to 300.0 sec. The second input on this card 
is the conputational tine step, 0.5 sec. The third input is the 
precision, 1.0. The second data card contains the initial conditions: 

• 

$ 0 , <j> Q and their scale values. The remaining data cards contain An^, An^, 
their scale values and the modal frequency, in that order for each mode 
being considered. Sample data cards are shown in Fig. 16 and input data, 
for a specific case, is shown on page 12 of the program listing (Appendix 
V - 12). 

The Flexbeam Output Program prints and/or plots any and/or all of 
the following: pitch angle (and pitch rate), modal amplitudes (and 

rates), and deflections taken at particular points along the undeformed 
beam, all as functions of time. 


(The designated points at which the deflections are to be calculated 
are specified according to the position of the particular point from 
the left end of the beam, non-dimens ionalized by the total beam length.) 

The desired outputs can be implemented by placing the "amplitude vs. time" 
deck of cards, (page V - 17 of the program listing), "print out-put" 
deck of cards, (page V - 16) of the program listing), or the "call 
to deflection plot (DFPLO)" subroutine, (page V -13 of the program 
listing) , into the main part of this program which is also on page 
V - 13. 


rUJJ iljw . • . *na ■ 
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6 . 2 Numerical results of Flexbeam Program 

The results obtained from the computer simulation of the dynamics 
of a flexible beam in orbit using the Flexbeam Program will be considered. 
The beam is assumed to be a long slender hollow tube made of wrought 
aluminum (2014T6). The length of the beam is taken as 100 meters * 
its outside diameter and thickness are: 0.05 meters and 0.005 meters, 

3 

respectively. The structural rigidity (El) of the beam is 7.707 x 10 

2 4 

nt - m and the mass per unit length is 9.906 x 10 kg/m. It is further 

assumed that the e.m. of the system follows a 250 n.mi (463.31 km) alti- 
tude circular orbit. 

The initial conditions which remained constant throughout all 
simulation runs are: d>(0) = 0.2 rad., {(0) = 0.0 rad. /sec., A n (0) = 

0. 5 meters (maximum value) and A^(0) = 0.0 meters/sec. Other parameters, 

1. e. the time interval of the numerical simulation and the generic 
modal frequencies, were varied. Variations in these parameters and 
their effect upon the deflection of the beam throughout its entire length 
will be discussed. 

The first two cases of the simulation were constructed to test the 
program. The first, depicted in Fig. 17, is a simulation of pitch motion 
with a superposition of the first and second generic modes, (refer to 
Table 2). Each node is assumed to have an initial amplitude coeffi- 
cient of 0.5 meters. The deflection is calculated using Eq. (11), where 

f V / 4 A 

(Fq) has only one component, k. 


The deflected beam was examined at three locations, the left and 
right nodal points of the first mode and the central nodal point of 
the second mode, (Table 2). In Figs. 17a fmd c, we observe the pre- 
dominant effect of the frequency response of the second mode, whereas 
in Fig. 17b, that of the first mode. Discrepancies between the responses 
shown in Fig, 17 and purely simple harmonic motion at a particular modal 
frequency, are attributed to: (1) snail numerical errors in calculating 

the exact location of the nodal points and (2) nonlinearities associated 
with pitch-rate coupling in the generic modal equations. In the second 
case, Fig. 18, the frequency of the first mode is set equal to the 
orbital frequency. The frequencies of modes two and three are calculated 
based on criterion consistent with free- free beams, (Table 2). Figs. 18a 
and b. depict the responses of the first generic mode and then all three 
generic modesj respectively. In each figure, we note the growth of the 
amplitude of the first mode due to orbital resonance as simulated. 

Fig. 19 is associated with Fig. 18b. This figure illustrates the de- 
flection of the entire beam for a one hundred second time interval. 

This interval was chosen to show the dominating effect of the first mode 
during this part of the response. It should be noted that the case shewn 
in Fig. 18 for the response of the first generic mode, duplicates the 
earlier result described in Fig. 5, after accounting for the nondimen- 
sionali 2 ation of and noting that the initial conditions in <{> and i 
shown in Fig. 18a, will result in a pitch amplitude of 0.2 radians. 



The third and remaining cases of the simulation shown in Figs. 

20 - 29, all depict the following plots; nodal amplitude and deflection 
vs. non-dimensional beam length as functions of time. The modal ampli- 
tude responses are used to study the interaction of the modes super- 
imposed upon one another; first with initial values of the node ampli- 
tudes equally weighted, i.e. An(0) = 0.S meters (Fig. 20), then unequally 
weighted, i.e. Figs. 22, 25a, and 27. Figs. 22 and 2Sa can be compared 
since the initial conditions are similar. It can be observed that with 
the frequencies chosen for the first generic mode, = 0.0628 in Fig. 22 
and - 0.628 in Fig. 2Sa. Which represents VlQO and VlO cycles per 
second, respectively, the responses shown, illustrate the effect of the 
greater rigidity in the latter case. The deflection plots are used to 
study the deflection of the entire beam during specific time intervals. 

As a particular example, we consider Figs. 20 and 21 jointly, to study 
the effect of the superimposition of the different modes at different 
intervals. We note that at the beginning of the simulation the signs 
of all three modal amplitudes (A^ - A^) are positive. Figures 21a. and b* 
show the deflection of the beam throughout the length of the entire 
beam during the first 21 seconds. From Table 2, it is observed that 
the deflections at the left end of the beam due to the first three 


(or more) generic modes , will be additive providing the modal ampli- 
tude factors have the same sign. This phenomena is apparent in Fig. 21a, 
where the initial larger deflection at the left end of the beam should 


be noted. 


The first and second modal amplitudes have negative values between 42.5 
and 48 secs, whereas the third modal amplitude has a positive value, 
(Fig. 20). It should also be noted that "die first modal amplitude 
reaches a maximum negative value at 45 secs. During this time period, 
the dominating influence of the first mode in the deflection response 
(Fig. 21c), is apparent. The contributions of the second and third mode 
tend to compensate each other. 

The remining amplitude and deflection responses showing the 
effect of different initial conditions and the numbers of inodes in the 
nodel can be examined in a similar manner and will be useful in the 
forthcoming simulation of the free-free beam under the action of 
various control devices. 

Figure 30 shows a typical response of the pitch motion of 
the beam for a simulation which involved only one generic mode. Since 
the pitch motion is not coupled to first order with the generic modal 
motion [Eqns. (88) and (90)1, this type of response is representative 
of all pitch motions simulated for small pitch amplitudes . 

For all numerical cases reported here, an average of 10-12 minutes 
computational time was '^equir* d to simulate the dynamics over an 80 
minute interval of real time using the NOVA 840 computer system. 





CONCLUSIONS 


This report presents the development of the equations of motion 

of an arbitrary flexible body in orbit. In the case of planar motion 

of a long, slender beam in a circular orbit, undergoing small pitch 

oscillations and flexural vibrations, the pitch motion completely 

decouples fit® the elastic motion and the elastic motion is governed by 

o 

a Hill's 3 - term equation. For large values of the elastic 

motion completely decouples from the pitch motion and the elastic notion 
closely approximates that of an harmonic oscillator. However, the 
numerical results indicate the possibilities of instabilities at very 
lew values of . A general computer program which treats both 

the pitching motion and the first-order flexural vibrations of a thin beam 
in orbit is developed; this program can be modified to simulate the effects 
of both external environmental torques and control torques that may be 
provided by actuators located at various positions along the beam. 


CO 
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Appendix I 


Transformation Relations : 
Transformation from to t 1 



Fig. (A-l): Inertial and intrinsic frames 

From Fig. (A-l) it is evident that. 


1 

H* 

H 1 


snow 

s n s u) c n 


r~ ~1 

X 


= 

c n c u 

c n s u s n 


Y 




ADDendix - II 



The transformation matrix connecting CP) and 7,, where P 
is a general point in the body, is developed as below, 

(i) Transform a vector in T-, CP) to a vector in Tq using the 
transformation ^ (?) . 

(ii) Transform the vector in Tq to a vector in t^CO) by the trans- 
formation (0) 

we car. write, 

T 1 ” 1 (P) = T 1 _1 (0) + AT n -1 

3Tn 3T-| 

= T. (0) + — An + — — Aw + H.O.T. 

1 3n 3w 

Hence , 

T.. (0) T^ 1 (P) = I + A (II - 1) 

where , 3T-, 3T-, 

A = T 2 (0) l-± An + ^ Aw] 

and from Eqn. (I - 1) 



Thus, the transformation ( I+A ) transforms a vector in t ] (P) to a 
vector in T. (0), i.e. 



II- 1 
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The gravity force per unit mass at 0, expressed in the 
intrinsic frame at that point, (0) is given by, 


f 0 (p, n,«) = vv| Q = 


1 9v 
psn Bto 


Where, V is the gravitational potential given by Eqn (5) . 

The gravity force at P referred to the intrinsic frame at P can 


be written as: 


(p, n, w) + -^7 < 

p 0 5 3(p,n,w. 


Ap 

An 

Aw 


+ H.O.T 


Hence, the gravity force at P can be referred to the intrinsic frame 
at 0, t ^ (0), through the use of transformation Eqn. (II-3), by: 

F = (I ' A) 7 

i,e * * ~ 9F r Ap ' 

0 3(p,n,w) ^ (XI -6: 

To a firsx order approximation in Ap, An, Aw, 


F = F- + A F n + 


o 3(p,n,w) 


FromEqns. vll - 2) £ (II - 3), it is easy to show that, 


£ F o 


-f sn 
w 

- F w c ’ 1 

F o sn+ F n cn 


II - 2 


where , 


F = K = _ va * ' + 
P 3p p2 


2 • va I K„ (|) S+1 


03 

' K 

S p 


S+1 


f = 1 jV - va ft3 . 

‘o p 9n p S 3 1 K s ( p } fi s 


F = _L_ 2Y. _ va » K -a, 
u> psn 3u psn - ^ ^ 


s=l 


. S+1 

a. ' 

S V ”S 


By some matrix manipulations, it can be shown that, 
9F 


- F 0 + 97p,n,to) 


Ap 

An 

Ato 


= B* 


where. 


B* = 


f£ I 1 / 3F P 

3p 


- F 


Ap 
pAn 
L psnAw J 


i / 3F 

- F 


p y 9h ny | psn ^ 9w 

3p i / \ 1 


OJ 


sn 


(II 


(ii 


9F 


1 ('* 


n j l l n + F 
9p p an p 


, /9F 

JL_ [ n - f on 

psn l 9w w 


9F 


w I 1 f Oi 

3p ‘ (> l 9n 


psn (9^ + F p sn + F n cn 


(II 


Evaluating the partial derivatives in B*, using Eqn. (II-9), we arrive 
an the rollowing ex ression. 

,2 


3- = 


va‘ 


) 


[S' ' + I K (2 : f B Cs) ] 


where. 


.( 0 ) 


2 

0 

0 


s=l 

0 0 

-1 0 

0 -1 


s p 


(II - 


- 9) 


- 10 ) 


- 11 ) 


12) 


II- 3 


,(s) _ 


(s+1) (s+2)fl„ 


-Cs+2)fJ' 

s 


-(s+2)§ 

STl 


-(s+2)J?g 


[n„ - (s+i) n] 

S 


~ i 

/M 
V ST y 


-s(s + 2)^ 

( M' 

Vni 

2 

[ft 1 cot r,-S7 (s+l+ '1 

s s s .-r) 

(II - 13) 


Also, to a first order approximation . 


ps^Aw 


(T 0 T„) 1 r 

o l 


(II - 14) 


where r is the postion vector expressed in the body axes frame , t^. 

After substitution of Eqns. (II - 10) and (II - 14) into Eon. (II - 7), 
there results. 


F = r Q + B* (T 3 T^)" 1 r 
Reprojecting on the body fixed frame, we get 


r = Tq + M r* 


where. 


f = T T F 
*0 3 2 0 

M = T 3 T 2 3* (T 3 T^)' 1 

If we write, H = + I K g (^) 


then. 


H (C) ^ T 3 T 2 5 (0) (T 3 T 7 ) _1 


(II - 15) 


(II - 16) 


(II - 17) 


(II - 18) 


M (s) = T„ T„ E (s) (T, T 0 ) 1 ^4 


The matrix M u is given in Eqn. (6). 
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AoDendix - III 


This appendix presents -the expanded forms of the vector expressions 

in equations (28), (31), (32), (43), (44), (45), and (46). 

Let us introduce the following notations, 

u (n) r c o(n) , 

H aB «a*B ^ 


r (mn) _ r , (m) , (n) - 
L cB *B 0111 

I (n) = i / <J (n) dm 

° m vol “ 

00 

By the assumption q= E A n (t) $ (r~), we have 

n=l 

q = E A (t) $ (n) (r n ) 

, n o 

n=l 


(III - 1) 


and q = E A (t) ^ n) (r ) 
*i n 0 

n=l 


(III - 2) 


Consider Eqn . ( 28 ) , 
«(n) 


E Q' 
n=l 


= / [r n x q + 2 r n x (to x q) + r n x (to x q) 


vol 


+ q x (w x Pq) - (r^ • to) (u x q) - (q • co) (to x r Q )] dm 

E Q <n) = E f [A (r 0 x $* n) ) + 2 A (r Q x (w x $ Cn) )) 
n=l vol 

• • 

+ A {r n x (to x ^ n ^) +(<^ r ‘^ x (to x r n ) 
n u u 

- (r Q • w) (w x F (n) ) - ($ (n) - w) (w x r Q )}]dm (III - 3) 
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Usir« the irtethods of vector algebra, the above expression car. be 
expanded in the component form to obtain 

i n \ + 2 k n [< 5 ♦ 

♦ A n [2<H^ * .£>> * - < + - (K - ) + K - > “= 


- 2 Vz (H « ' ' Vsr tH - + H - ’ 

+ u u> flfi 5 * *&>> + <«- 2 - V> ( <z’ + H z? ))] 

x s xy y* * - 


cm - 


0 (n) ^ Q <n) ^ obtained by the cyclic permutation of x, y, z 
z 


in the expression for Q 


(n) 


Now consider equation (31), which is given by, 


Gr = / r Q x M r Q dm 

vol 


On expanding in the component fora, 


(III - 


Sr = (J z - Jy) M 23 i + ( J x - J z ) M 31 j + (J y ” J *> M 12 k 

where, M. . is an element of M matrix 

Jx , Jy , £ are the principal nvaments of inertia of the body in the -de- 
formed state. 

r» 1 


Considering Eqn. (32) for G n ' , we have, 
= / (r n x Mq + q x Kr Q )dm 


® (n) 


I G 
n=l 


vol 


1a/ (r 0 X M f ln) + H-> (ri) x M r 0 )dm 

n=l n vol 


(III 


III- 2 


Hence, the components of (III - 6) are, 

G x n> = V m 33 - «22> (H yS 5 + H zy^ - « 2 1 ^ + H zx >5 

+ M 3X< 5 + ^ + 2 M 23 «$’ - O 
G^ and components are obtained by the cyclic permutation 
of x, y, 2 in (III - 7). 

Now consider the following scalar quantities. 

Eqn. (43) : 

4> = / [$ r '^ • to x r n + • to x (to x r n )]dm 

n vol u u 


T (n) TT (n) 


= “x G£" - H-n ♦ o, y (H-' - H-') ♦ i z ($> - $>> 


+ Vy + ^ + Vz + H zy >5 + “A <H « + 

- V + H zz }) - V «zz’ + H ™> - “z 2 «£’ + “w 5 (III - 8) 


Eqn. (44) 


_ , 'inn , 
m=l vol 


= / [2 • co x q + a x q + • to x (to x a)]dm 


= I / [2 A (* (n) • w x * (m) ) 

m=l vol 


+ A^(<i^“^ * w x • to x (to x $^ m ^)]dn 


Hence , 


mn 


nr 


(mn) 

J 

T (mn), , n (mn) 

- L ) + to (L 

T (mn), 

- L ) + to 

a te) 

- L (mn) ) 

yz 

zy y zx 

xz z 

xy 

yx 

(mn) 

t Cm), , * n (mn) 
- L . ) + to (L 

- L ) + to 

( (mn) 

- L (mn) ) 

yz 

zy y zx 

xz z 

xy 

y x • 


n (mn) , T (mn), , , T (mn) , .(jnn). , , T (mn), T (mn, 

+ to to (L + L )+toto(Li +b ) + to to (L + L ) 
yx yz yz zy zx zx xz 


x y xy 
. (mr 

^yy 


- “x 2 + - “v 2 <?' + £”> - “z 2 ^ + 

(III - 9) 


III- 3 


(Ill - 10) 


Ecn. (45): g = / ^ • M r n dm = II M « 

^ vol 0 a B aS aB 


*-(n) 


Eqn. (45): I g^ = / $ • M q dm = I A^ J 

m=l ‘ vol m=l vol 

Hence, 


r $ (n) • M $ ClIl) dm 


g. = A I I M 

°mn m « ag ag 
a p 


where, a g = x, y, z or 1, 2, 3. For example, when a is x in , 
the corresponding value of a in M^g is 1. In a similar way when a is 
y in H^, a is 2 in and when a is z in , a is 3 in Same 

reasoning holds for 6 also. 


(Ill - 11) 
,(n) 
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.•'Appendix IV 

Natural frequencies, mode shapes and nodal mass for free- free uniform, 
beans : 

The natural frequencies of a free-free bean can be obtained by 

9 

solving the following frequency equation. 

cos Si cosh B£ = 1 (IV - 1) 


where , 


i - length of the bean 


s 4 = 


arm' 

El 


m’ - rass per unit length of the bean 
El - bending modulus 
a) - natural frequency 

The mode shapes corresponding to the frequencies obtained by Eqn. 
(IV - Dare given by, 

[/ cos& n £ - cosh S n £ 


Z ,(n) (xD = D 

n 


u 


sinhS i - sin S i 
n n 


(sin 6 x* + sinh 6 x 1 ) 
n n 


+ (cos 6 x’ + cosh S x’) 
n n 


(IV - 2) 


x 


Fig. (D-l): Free- free bean 


IV-1 


where, g_ (n = 1,2, - — ) are roots of Ion. (IV - 1). 

It should be noted the” is measured from one end of the beam, 

. £ 
as shown in the Fig (D-l), Hence, x’ = x + /2 

Define, 


- = c 
£ ^ 


(IV - 3) 

Hence, the non-dimens ionalized form of the mode shape in Ea. (IV -2) 


S£ = £2 
n n 


is. 


Z Cn) ( o 


cosfi - coshft 
n n 

sirihfi - sinft 
n n 


(sinfi^ £ +si hhfi 5 ) 


+ (cosfl n £ + coshft c ) 


(IV - 4) 


The following Table - 1 gives the first five natural modes and the 

approximate values of the corresponding natural frequencies of a free-free 

. 10 
beam. 

The modal mass (generalized mass) in the n^ mode is given bv, 

v2 


M = / 
n v 2 


(x) y dv 

Since the beam is assumed to be uniform, y = constant. 


(IV - 5) 


Hence, + £ /2 

M = m’ / (x) dx 

x= - £/2 2 


CIV - S) 


Using the change cf variable ; = + y , we get 

M = m'£ f [I (n) (0] 2 d; 
n 0 

Substituting (IV - 4) in (IV - 7), 


M = D " m'£ f [K (sinD £ + sinhft O 
n n 0 n n n 


(IV - 7) 


+ (cosft C + coshft OKdC 
n n 


(IV - 8) 
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Table - 2: Mode shapes and frequencies cf free-free beams 


Mode shape 

l-*- c 5 = 1 


rrequency 



oj 1 = 3 . 56 fZl 

Vn'r 


\ 

\ 


0.13 


\ 




— ~ 0.868 




\ 

\ 



\ 0 . 09 ^ 

-V— 


0.356 



0.644 0.90^/ 




^0.0734 .277 .5 



.723 .927 

\ 



\ 0.06 .227 .409 .591 .776 .9 





= 47 



Note: The frequency value of zero corresponds to the rigid body 

nodes. 
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where K = (cos ft. - cosh ft )/(sihh ft - sin ft ) 
n n n n n 

The various definite integrals appearing in Eqn. (II 

evaluated in Table - 3. : 


- 8 ) 


Table - 3: Useful definite integrals 


r a 2 n »- A* - - sin 2ft n 
f sm ft C at; - - - — _ 

0 n ^ 4£^ 


f sin ft C cos ft £ d£ = s '*~ rl ffn 

0 2ft-, 


/sin R s Sirih 0 C dc = cosh % sin 0^ - sirih cos 
„ n n 


2^n 


/ sin ft r cosh ft £ dr = s:i:nl1 ^ “ cosh cos ^ - 

n n n ' “ 


2ft. 


n 


/ cos 2 ft C d5 = i + SJJl 2S V 
0 n 2 4 


/cos sirih iy; d5 = cosh ^ cos 1^ + sirih fl n sin Sfr-1 


2ft 


n 


/cos ft c cosh ft ? dC = ‘ sirih ^ COB + gosh fti sin ^ 
n 


2ftL 


n 


/sirih 2 ft C d? = 5ijl h 2 fo 

ft n 


4ft. 


1 

2 


n 



are 


1 

2ftn 


/ sirih Q C cosh 0 C d£ - sirihO /2E2 
0 n n n n 


/ oosh 2 B s d C = sijh 2g n + I 
0 n 2 . 


40 


n 


Substituting the values of various definite integrals in Eqn. (IV - 8)? 

after simplification we get 
2 

v — ^ p/i*- 2 ■ ^ sirih 20ji . ~ » — . 

K n - 2 L [ 0 ^ + 1 ) (-— — a + cosh a n sin ) 


0 


n 


- (K^ - 1) ( S .f n * sirih 0 n cos & n ) 
+ K n (sin + sirih + Oj 


(IV - 9) 


IV- 5 


mu ^ffWCSiabww*-* 


B 


HDWARn UN 1 U£R-S I T Y 


SCHOOL of ENGTMFFPING 


MOV ft Rao computer 


* 

* 


£ 


*********** ******-***************** ****** ***************************** ********** 


PTRAN/P/E/P T m p TfJLRR/nMP.LS/L 

S COMP I L £R 'DOUBLE' ’PRECIST ON 

FLF Y PEft^ main PROGRAM 


REAL Lljrin, L31 rL 33 f m u y 1 3_rj*31 ,± 133 ,MS, J2 ,,J3rL SjLJLNU. 
’COMMON /BEAM/ PARMCn / Y(02 ) ,DY ( 42 ) , 

2 HI lC20),H33C20),Hi3C20Y, 


I _3 L \ 1 00/20) ,L3IL20,?O^Lli_C?0 ,?0) tL^lf.PTLrZPAi 

' , ' i- ii r.MSl , Ms , ALPHA.* 

5 PC20) ,PPC2Qr?Q) *GC2<n rGGC 20-,20 ) , 

6 _DWf 20) /H.RC20) r HDC?Oj , 7 , MU , A , RILJ2 , J3 , 

7 0CC?.Q)/6cP r .«6/EC2d)rP2fLS'C20f 20.V#LDC20f 2*6) ,C2 ,m,nqRD 

COMMON /LOG/ SIZN0MC4P) , SI7MAXt«2), LIMES * LSTEP 
DIMENSION ft U X C 1 2 / tl2) 


1 OPEN 5,”$YS0UT", ATT="AP‘' 

* 

i DELETE “FLEXSEAMm 



OPEN 3, "FLEXSEftM'f 

f GFE N 2, 11 I MP 11 , ATT ^”! ” : 

CALL FRINP(PRCRN) 

CALL FBSETCPRCSN) 

B LIN ES = 0 

LSTEP “ 0„20*C ■ CPARMC2T - PARM Cl ) )/PftRMf 3) ) + O.S 
* s CALL PK-GSCPA D M/Y/DY r NORD, IHLF/AUX) 

DO 7 1 = l r NOPD 

Ij RATIO = CSIZMAXCn/SIZNOMCin^lOOi.O 

P WRITE fS/ 91) tr PATIO 

7 CONTINUE _ 

H WRITE C5/9) IRLF 
11 CALL EXIT 


U 9 FORMAT! r OFINAL VALUE OF THLF ; ' , T 3 ) 

91 FORMftT_C f OYCJ^r 12 / ') PEAK ; * ? Ffe„n 
„ 'end 
tIplist 

-pr[NT 1 _ _ __ 
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FFJLETE TWP.LS 

criTPtM/B/E/P T W F TMP.LS/L FCT.RB/* 

® COMPILER DOUBLE PRECISION 
SURF OUT INF FCT*' CT, AY^OYl 

I REAL LlULn.L^ULSR.Mt t f J?,J3,LS,LD P WU 

« COMMON /BEAM/ PARMfST, YC«?IfOTCa?) , 

? "Hl'i ('SOUH^SCpOVrHlJfaOTf 

S 3 L 11 C?Cl,?m , L3M?0,20) , L 1 3 C20 , 201 r 13 i C?0,20>, 

6 Mi 1 »M33,Mf 3,M31 ,MS, ALPHA, ^ _ 

S PC 2 6 1 7 P P ( 2 6 /go") t G C ? 0 }7 GG C 2 O , 2 6 ) » 

1 6 OlvCcO) ,HSC20) ,HD(20) ,?rWU, A,Bn,J2,J3, 

7 QC C2Q) ,QCP r WO,E C201 , R2 , LS C 20 , 20 1, L D C20 , 20 1 £ C2 f_M,NORD 


i 


FRPIF — DTFF EON'S 


CALL FRPIF 
CALL FBSIF 
DY C 1 5 =ALPHA 


rjYC?i=Ycn 

DO ? 1=1, « 

_ DYC?*I + n=FBAPRtI) 
DYtP*f*2)=Y (2*1+1) 
CONTINUE 

RETURN _ 





f VPLIS'T 
F?IMT 1 
RRT^T J 

jflPTP'A^/a/E/P T‘.‘P rwp-LS/l FRRET r RR/R 
1 COMPXUxP DOUBLE PRECISION 
SURROtiTTME PRSETf Q PPC5WV 



I 



PEAL LXUL n,LM r L 3 r M t 1 , * i ^ , M3 1 , , ms f Jp , J3 , LS , LD , MU 

COMMON /ream/ p A R m ( El f v ( a 2 ) ff D Y C A ? ) , 

? H 1 i C ? 0 ) f H 3 3 C ? 0 ) / H t ^ f p 0 1 r 

3 LH Can,?OD f L3 vCao,2nWL1 3C?0r?o) ,L3t C?0, ?0), 

« ^11 fM3-5 r Mt-3 r M^ i,ms, ALPHA, 

5 PC20),PPC?0 # ?ft)rGCPO),GGC?n,?Ol, 
h .Divf.203 *H$C20) ,HDC? 0 ) ,7 r MUr A,RO, J2» J3, 

7 nccacn ,RCP tf «D,Efan),R? r LBfan f ?n) ,Lorao,?o) ,c?yM f nord 




PBS FT "SET ‘IP TMTTAL STATE 

CPF?ECSN=PRECTSTPN, RELATIVE TO t.O) 


’CALL FBHST 
CALL PfiLST 
2 sI. 2 «iUE -6 _ __ 

CALL F>SCL CN6rD,OY # OV,p*RECS'n,PAr‘m) 

return 



END 







ORIGINAL PAGE IS 
OF POOR QUALITY 





agteaaMgffiateiiteiaaia^aaBiifflfittBffiaga^ — 


jUttBasa <»■-- «-*>**» 


fifim . 

| NT \ ' ■ ■ 

,T»P T«P,U/l, FAHPT.R9/R 
CnMPUFR POUBLF P*>FDiStO w 
SI IMRQUT I M® tRHST 

PF At Lit pi i 3, L ? 1 » L >3 , MV V , M rs f M3! t , i WR , J2 r J 1 * t.U> l&> WU 
COM^DW /BHAM/ P&IJNACS'I r V C ^ ? 1 
5 Hi 1 f?{O,«33C2C0 ,HHt?0V ( 

3 U..H?0,?0),LX^(.20r?.n^L.1 , 5f»0.,20WL\t {.?n,?03, , . 

<1 Ml l f 1.'^ f M13 v ,MS, ALPHA , 

5 Pt20>>PPC?0r20) ; »Q(?ni / Gt?C?n-,?0)» 

6 pi*f ?0) ,WSf?0) ,HDtPO) rZpNMI, A,Rtl, J2, J3, 

7 QCC2M pQCPp .‘in, FfPAt ,RP,i?(Sn f ^o'i pLnC?r«,? 0 ) f C 2 pW,MORD 

FRHST — SET *HTO r 

DO. 10 I-1,M 

HI l,C I !) ” 0 * <■*■- . . ...,.' ...... — la.'. ... a...— 

Hismso.o 

H33CI)pQ*0 

hs m=o, __ - ~ 

HI) (1 5=0.0 

to cnwriHUE 



to 


PMGSNSE PAGE IS 
OP POOR QUALITY 



t-a/» 


wptxsr 

■PINT 1 


IFORTPAM/p^H/p TMP TWP.USVt FRL5T.PR/R 

S COMPILER noil RLE PRECISION 
SUBPnutIME FRIST 

REAL Lit ,Ln f llt ,L33r*Ml r M nrM5t,M33;^MS 
COMMON /ream/ PAPVCP) t rOYCaPT f 


, J?r J 3 *LSfU 3 rNU 
?V 20 3 , 


I 


? Hi t ( ? 0 3 , H 3 3 C ? 0 3 r M 1 3 C 2 A 3 , 

3 LI 1 <?0,?0) r L33(?n,?0 3,U3C?OrH0 3 ,L31C 

a MI 1 ,M 33 ,M 15 ,M 31 ,MS, alpha, 

5 PC 203 r PP(PO, ?n) ,Gt?nt,GG( 2 Cl f ?ni , 

6 DWC 2'0) fHSCPOt , HD (20 3 , 7rNtJ, A,RQ,J?,J3 r 

7 DC (20) ,QCP,*D,Ef203 ,R?,LS (20,203 ,LDC20,?0) ,C2 ,m f nord 

FBLST--SET ' LT J ' . 

oo ho r=t,M 

DO 20 J=i ,M 
io in ttrJ)-n.6 
00 30 1=1 f M 
DO 30 J=i,M 
30 Li3(irin=o/o 


«0 


61 


DO 

DO 


40 

AO 


I»l- ,M 
J= UJtL 


LJt C X r J)=rft„0 
00 50 . 

DO 60 J=t.,w __ __ 

L33U,j)=0.cT 
L 3 S'Clf I ) “1 • 0 

DO fel_I = l,M 

on 6 i j~i r m 

LOCI , J3=L3t. (I t J)»L13CI r JV 
1 3 ( I , J 3 =L 1 1 f it J 3 +L 33 C I , J ) 
continue 

RETURN 


END 






LIST 
NT 5 

RTRAN/R/F /p TMp TMP,LS/L FRFIF„RK/R 

[ COMPILER r>nURLE PRECISION 
surhouttnf fb d if 

p F & L L It , L 1 3 f L 3 .1 , U Wit r M 1 3 1 , MJ 3 r y P , J?> J 3> LF , L 0 r W1J 
COMMON /RgiM/ PAWR), Y(4?) ; OYt/l?) , 

? HI 3 C205 f H33(?0) rHISCSO), 

1 3 L 1 1 ( ?0 f 20) * L33 C20 j« go I f Li 300 ^20 V r L3 1 C20 r?G ) r 
u Vil/M33rl*U,MM#MS*&LPWA» 
i 5 p(?0) f PP(20r20) P G(?0) ,GGC20,?0) *- 

\ b Dw C 20) ,.HS.C20)rHDC20Yf7rMU f AirROf JP r J3» 

7 GCC20) VQC p , ‘iD,F(?0),R? r LSC20, 20) ,LDC?(V,?0> ,C?,M,-NC]RD 

frpif-^primary intermediate functions , ; 

SPH=STN CY C2) ) 

CPH=C0S C.Y.C?) ) 

Mil =3* C PH *■* 2 '-V 
M13=-3*SPH 
M3U-3*SPH*CPH 
‘ M33=3*SPH**g-'r 

MS-MI3+M31 
p?=vcn+‘wn 
RETURN " 



cMn-UA 




•fja \ 


f MPLIST 
U W t 

IFOPTRAN/P/F./P TVP T^ c .LS/l. F«STF * RR/B 
M COMPJLFS pnu*L£ PRECISION 

1 SUBROUTINE FRSTF 

INTEGFR « 

$ RF-*L tl ULnrUl,m,M1 J?,.I3.LS,L0^lt 

£ • • C0MWC1N / B F A f V PARMCSUYfU?) r - ■ 

3 H1 1 C?0) r 

1 3 m c?a, 20 ) # L33C?o,?m ,Lt3r2o,?n)>L3i r?6, ?m , 

« ■ ^1 J fM33rMj-5,M31 r WS*6L.PHA> 

5 PC?o) f epf?o,2ti)»«c?ov,Rt;c.?n f ?fty, 
b DLC20) ,W$C20) ,KOC?n) f 7i MU, A,RQ, J2, J3, 

1 7 .QCf?O)*.QCP>rflDi k Ef-?Ol r Pa,L;-S-C?0*?n)fLpr2Of?n),C?,^,WOBb 

C FRSIr— secriNBARY TMTFRMpn T ATF FUNCTIONS 


Y2~Y (2) : , ' _ ■ 

”TY2=2«0*Y2 ' ~ ~ :'“~™ — 

CTY2-C0S (TY2) 

$7Y?=SX^CTY?) ■ ■' ■ 7 - / 

SUM1 0»0 ~ “™ ■’■* 

5UY-2 = 0*0 

_ „ . DP , 3 R “ 1 , M • .. 

SUWi = SU^l 4 3 /n * f H 1 3 f*3 * GT Y ? V" & » 5 *»%( R ) *$ T Y ?. i * Y C 2 *R +? ) 

* - HSCR)+R2*YC?*R-M ) 

$U M2 = SU*2 + HS(R1*YC2*R'*S1 

" CONTINUE ' V ■ ~ ~ ~~ 

■ALPHA- = Z*C-1 .5*J3*STY2 +'2.0* : SUM1 ~ C2 - GCP) / CJ2 + 2,n*SUM?) 
DO 70 R = 1 ,M 

PCR3 i=-R2 : fc*2*HSCR) ’ ~ 

G(R)=Z*CHU rR)*«U+H?lfP>*MTSl+?*HtlCR)*M513 
DO 80 Wsj>«- 

PPCR,N) = (?*YC2*M-H)*R2ma*N+2)*ALPHA)VLDtR,M) 

2~Y C2*N+2) A-R?*ft?*l.S(R f N3 

GG C R , N 3 = Z * Y C 2 * N +? 3 * ( L 1 U R , N ) *M 1 1 + L 1 3 CR , M ) *w 1 3 +L 31 f R , N 3 *M3 1 _ 

24L33fR r N)+Y33) 

so CONTINUE 
70 CONTINUE 

RETURN ' * • ' "" •■• 



<PITST 
fl WT I 

JPTPAN/P/f/P TMP TMP.LS/L F**A2R.RB/R- 
■ COMPILE*? DOUBLE PRECISION 
FUNCTION FBA2=JC*?I 
. INTEGER R 

PFAl m ,L1 3f L3t ,L33,'-M 1 , Ml 3 , M3 5 , m*53 , ms , J? , J3 , LS , LD , MU 

COMMON /BEAM/ PAPMCS) ,V C<I2)*DYC42) , 

? Hli(2(n,H33(2GI,Hr5f?<n, _ 

3 LI i'(?0 > 2 6 ) f L33 C2 0^ 2 0 )>Ui 3 C20 / 20') j»L3'1 C20 p'20 3 / ■ 
fl MJ J >M33> M i3,Mli'. r MS, &LPHA, 

5 .P£20-),PPCPO f ?n),rj(2ny,GGC2'n,?0"5r 

6 DWC?0) r HSC205 r HDC?flK?fNU- r A,RQ, J2, J3, 

7 NCf2n/, QCP, wD,E-C20i r«2f LSC20V20') .LDC20 r ?0) ,C2' r M,M.0RD 

FBAPR--EVALUATE '2«ND DER OF A(R) " 

SUMnO.O 

bo 

SUMsSUM-pp fp,N'HGSfRfW) 

Qn CONTINUE .. _ „ 

F8A s PR = -0WfR f**?*Y.(2*R+2) *P C R ) +SIJM+S CRT +E CP I +QC (R ) 

RETURN 


ORIGINAL PAGE IS 
OETOOR QUAMT’ 


p 



tg^UST 

jj|»lNT 1 

ILF? N DP();FUFXPEft^5*LH/n 

* m OUTP.RB .FBSET.R* FPHST.RR FRLST.PB FRPTF.PB FftSJF.RB FBA?R„Rfi 
ilELETF 

OUTP.PR FRSET.PS FBHST.RR FBLST.RB FBPTF.RB FRSIF.RB FBA2R.PB 
IGTOD 

llELETE TWP.LS 

iIqpTR A<M/»/E/ p T^P TMP.LS/L OiITP.RR/8 
COMPILER DOUBLE PRECISION 
| SUBROUTINE OUTP CT , AY r AOY , Th*LF ,N, AP ARM) 

X PFAL L 1 1 , L 1 3 , L 3 l , L 33 , M 1 1 , M\3,M3i ,M33,MS, J2, J3,L5 r LP f 'HU 
COMMON /RFAM/ PARMC5’) , YCt|?5 r DY(/4P) , 

f p Hll can 5 r H33C?0) ,M\3f?0) , 

3 111 (20,20) ,L33C?D,20) I-L13C20, 20) ,131(20,20) , 

ti MH/M33,Mi3,W3 5 ,m$, ALPHAS 
5 PC20) ,PP (20,20) ,GC?0) ,GG( 2 0,20) , 

_6 DW c ? 0 ) , HJSC JJO)/ HD C2 0 ) , 7 . , MU, A ,_R0 , J2, J3, 

7 DC(205 ,OCP,WO,EC2Q) , R2,*LS C?0,?b ) , LD C2"6 , 20 T, C2 f ~MOPD ~ 

COMMON /LOG/ SIZM0Mf.«21> SI7MAX («2) , LINES, LSTFP 
LOGICAL pknyt 


r 


F80UT --TAKE OUTPUT VALUES 

"OF 2" ”l~" i V i”7~~n6rD “ " ' ~ 

YM AG = ABSCY(I)} 

IFCYMAG ,GT, SI7.M4XCI) ) STZMAX(I) a Y^AG 
2 CONTINUE ' 

IF (.NOT. RKNXTCIHLF) ) GO TO 8 

WRITE BINARY (3) T, THLPjJTY Cl ) , 1=1 ,_NORO) 

IT = (T T - PARM'cn 5 ./PAPMC2)i*io6ro + 6.5'“ 
TF (MDDCLINFS,LSTEP) .'ED. 0 ) TYPF t't,THLF 
LINES - LINES + i 
8 CONTINUE 
RETURN 


% 

. 1 

Ri' 


: TRAN/R/r/P T^P Tn3 r RR/R TMF.LS/L 
COMPILER DOUBLE PRECISION 
SUBROUTINE PPTWP CPRCSN) 

PEAL 1 1 1 ,U3 ,l\X rl>*r «t 1 r «l 1 r ***' 33 * *S .. J? , J3 r LS, U* * m 
COMMON* /SEAM/ PARMf^ ,y ftfpunYfilRy, 

? pmpoXrHiXrsfiUBm^on 

? Li t C?R,20),L3K2a,?0),U3C?0,20),L3t C20>?0), 

0 MnrM*53# w t3rM^UM‘3 r &LPHA / 

5 pf?n) , PP(? 0 r 2 O)#Rt 2 O)fGGf 20 r? 0 Jf ■ . 

6 DWI20) ,HSC-20i »HO'C20> r A r ; RO, J2 f J3, 

7 RCC20) ,OC^p«'0fEC20l ,f*2>L3C?0 f 20) *LO(?Qn?03 ,C2 F M,K»nRn 

COMMON /LOG/ SI2 r *mMCapXr SI7MAXC02)r LINES, LSTEP . . 


FLEKBEAM 


imp i.rr 


MMP1 = 2*W + 1 
wypg — 2*w + 2 

READ f?VPl) CPARMCH, TjsJUJJU FRCS.M. ._ 

" WRITE "(Sr 91)'* (PARMCI), I = 1, 3), PRCSW 
READ (2 r *2) Y ( 2 ) , Yd), DYC2)* DYCU 

WRITE. f5r R2),_Y U 5 ) Yt I ) r_„DY (?) f . JDY t\Y ... 

' DO 2 L" =' \ t M 

II — 2*L F ? 

I? = ?*L +_J . . .. . ......... .... 

READ" (2, B2) YCTUr YCI23* DYCt!), H Y ( T? ) , OW CL 5 
WRITE CS,RX) L f YCI13 r Y-CI2)f D-YCTl)r DYCI23 * DWCL) 
2 CONTINUE _ . ... 

do n i = n moro * 

SI2MOMCI3 ~ DY Cl) 

SX2MAXCI3 » 0 o 0 ■ 

a CONTINUE 
RETURN 


Ibi FORMAT CAF10.Q) 

0? FPPWAT C7F10.0, FR^fl, IX). _ '_. 

<U FOPMAT c'OT FROM', F7.1, f TO', F7.W # RYS F7,l r 6X, 

I 2 ' PRECISION F<J„«0 

FORMAT ctnXf 'INIT VAL IN IT DER jSCALF UAL. SCALE. DER’, 

2 6X, # OMEGA'" / v THETA % ?FH,6, Ft 3, 6, Ft 2.6) 

1 93 FORMAT C* MOOES I?, 2FU.&f Ftl.&f Ft2»6, F13„6) 

i'TMPLI$f " 

; STINT t 



£ m -ft V/P./E /P r^R TO?. Rft/R T VP . LS/L 
CO^PUER nmmLP PRECISION 
FLOCK DATA 

REAL Lit ,U3,L3l rLSS^Hlr.^tSr^?! r v -T>, w S,J?, ,J3 r LS„LD,HI 
COMMON /REAM/ 1 , Y t ttg-J , DY W > , 

? HU f?f>) ,Hnt?0U 

3 Lt I f?n r ?0UL3RC?0.?f>ULt ^C?0,?nUL3t f?o r ?OU 

4 Mil ,v^ r M}X r M3t,MS, ALS'HA, 

5 PC?OU p PC?0^50UOC?Ol/OGC?Or?OU 

ft D‘.n?0)>HSC?03 ? H>)C?0} . f ? r MUr4 r PQ r i? r J3> 

7 RC(20) ,RCP/^r.i,EC?0) j.P^eLSfJO, ?0) »>LD(?fl,?nuC2f M *NCi p n 

COMVqm VRKMVlJ/ IHQLD, NE\ f T 


DATA P/E/, WbRD/lg/ 

DATA HU/9, 7F7/,A/637R'?R8./>R0/6FfUiBR./,vT?/6?55007/ 

data J3/85EH0 0 , / / QCP/O % O/f '40/, 00 1 1 15/r C?/0 , 0/ 

DATA Y/Q2*0» 0/ __ _ 

D A T A 6 Y / 42 * 0** 67 

data otf/ao*o*6/ 

DATA QC/0.0U9+0./)/_ ... ... 

DATA E/o'.O V 19*0,0/ 

DATA PARM/5*0 . 0/ 

DATA IHQLB/O/, HExr/1/ _ 


END 
«Pl 1ST* 
tNT 1 


-0KDB3METAGE IS 
OPPQOBQGAIOTf 


■ Wv > 


HOWARD UNIVERSITY — SCHORL OF ENGINEERING — NOVA *40 CDMPtjTFR 


\*&?** **************** *********** ***- k*ic *** * **** fr +*****^ r ^ ****.* ** 


| ^ FLEX RE AM INPUT DATA 

OP C<; FLEX REAM " " 


300,0 
0,0 
o.o 

}« j | o,o .o 

• 3 l . _ 0 . 0 . .3 

1 2 ^ 0.0 .? 

! t «■ 0.0 , 1 . 

GfH > 

jofLETE DPO ;FLEXBEAM 
! XFER FLEXREAM DPO ; FLEX BEAM 


1,0 

. 00030 

.03 

.003 

.003 

,003 

.003 


.06203165 
,1733170607 
, 3370 £ 7 _ 05 _ 
.5579 A 68 
,6337786 


\ 

T 

T 


~ — 0 £'-£ QOR '- 


..*.*•*••*>*♦ *:*..**.**:*.*•** * * *•*.-* *+•#.***.*;****;* *;* ft-fr+Jr* *■**■+** *■*■* IPtr* *••***•****•* **.*'*.***.•** ****** 


HOWARD UNIVERSITY — SCHOOL OP ENGINEER J NQ NOVA r/iq COMPUTER 


^ * > +• * >'*, * ;*•'* *,*;'* '*' * * * * A *** ***•***-*.*•*:*•* ****• *.* *•* **r ******** * * * * *•*'* ***** ** * * * * * * ** 


ir_ FLFJfftFftW OUTPUT PROGRAM 

//“ FOR TRAM ' , ’* 


1 PARA^FTFP MORO ” t2r POX MTS = 601, MOOES = 5 
■ DOUBLE PRECISION TT, YCMOROl 


DUUSLF PRECT S IJTN TT , YCMORIXT 
"DIMENSION ff POINTS)’” 
dtmfmsiqn PHI (POINTS) 

D IMF MS TOM AC POINTS, MODES) " 


JCALL IMOJJT • (?, S) 

CALL FDPEN (3, 11 FL EM REAM" ) ” " 

DO ? I “ 1 , POINTS 

PE AD BINARY C3V IT, IHLF , Y 



ACT, J) = SNGL CY C2*-*Jf2 j ) 

1 CONTINUE 

CO NTINUE ' 

CALL DFPLO (POINTS, MOOTS, A) 

CALL EXIT 


!<V FOPTP’At'J 

I ' WftXM « Sii MX- S. t o t 

“ PAP&MEtFR IT I a Q5> TT2 a H6, DT a a 
PARAMETER ftMAx a 3,0 
1 SU P R OU T I N E D F « L 0 C NPTS ,. MOOFS* AT 

I DrPFMSXa^ ACNPTS, MOOES), SHAPE fMY M4yw) rDFF(MV) 

DOUBLE PRECISION OMEGA 

f COMMON /nFPLDV OMFGA(MAYM) 

DATA V4, ThO nan ado r 7 , AS35PAAD0 , 10,993607^00, 

? ' 1«.n~t655U0, t7.?7E7S96D0/ 


1 


f 


a 

s 


CALL PFPLS CO^EGA, MODES, NY, SHAPF) 

CALL PSIZE C5 , 0 * VI .O') 

CALL PFRID Cl, in 
CALL PLDGD CO, Or 11 .01 
CALL PvSIZE CS,0, l-.OT 
DO e I T a m , TT5, DT 

1 ‘call plogo *co,o7 Cor 

on 5 IT S 1 f MX 

OaO , 0 

' DO a WOD%$”*‘ ~ 

& K Q - + ACTT, I^TASHAPE Ciy,TM5 

CONTINUE , , ; _____ 

OFF (TXT a Q 
CONTINUE 


: CALL PLOY t«DFUX> OFF, DMAX, NY 5 

W C ON T I N U E r rr~‘ ' • - t ~ ‘ ~ 

RETURN 


-two 



-1 


// FORTRAN 

f SURPOUTT’VE DEPLSCO M EGA, ^GOES, NX, SHAPE) 

DOUBLE PRECTS.TOW o^EGA C w OOES) 

DIMENSION SHAPE CMX, MODES) • 

§ DOUBLE PRECISION D7 , -W, EW, SHW, CHW, COEF 

DOUBLE PRECISION ZETA, W7, £W?, SHWZ, CHW?, DSHAPE 


= 1 .ODO/DFLOAT f NX - 1 ) 

6 .-.XM s 1, MODES 
W = n«E6A (XM) 

FE = DEXPCW) 

SHW = CE.' J ~ i * 0 D 0 / E W) / ? . 0 0 0 
CHw = cew + U0Dn/EW3 /J?.0Dn 
CP£F s COCOS CW) - CHW) /(SHW - OSJMfW)) 

. DO a ix = 1 > NX 

ZET A = DFLOATCIX - 1)*DZ 

WZ s, W*ZETA 

EW7 DEVPC'WZ) 

SHWZ = CE*Z ~ i.ODQ/EW7')/?.000 

CHWZ s CEWZ * l » ODn/E^Z ) /2,0OG_ 

" D S w APE ' e C 0 EF * C D S IN C W ? 3 + SH WZ 3 + ” 
SHAPE C IX# IM) = SNGLCDSHAPE) 
CONTINUE 


COCOS (W?) + CHWZ ) 


8 CONTINUE 
RETURN 


J ■ END 

1 LOAD PL0T./L 

f%Q dpoTflexbeam 
I XEO 


FLEXBEAM 



I 


I 


DO & L - i , HI 

I PH T =t 1 0 a-L « 9 / 

I At A = L ' A'^OPV-' 'v ■■■■- ■'■"■ ■■ ; * ' ' ’ 

I A 1 B = L + 9 ftQ 

4 CQMriKMtc 

DO 6 L * A?, lot v~ , v ~ 

1 PHI = 1 CV*L - 9 
ft CC^jTIHUE 

t* ; P I TE C5,9n TCIPKI1 p PHKTPHm' T'ciil A)', A tlAf* 1 } T A'C T A 1 aJT) 
X - * T C T A 1 B T> A Cl A j B r 1 ) f AfIAtR,i) . 

tv’PTTF f5,9H TflPHT), PHTCTPHI’j 
91 Format (\y; FfiVf, F j 4 * 9 P ( F 1 3. l', i?F 1*4 . 9)*) ' : 



ORIGINAL PAGE IS 
OP POOP. QUALITY 




PLOT AMPLTTUnE VS. T7^£ 

CALL PSIZE Ch P 0 P 5>0) 

CALL PRO* 

CALL PAXES 

CALL PLOX 00.5, ACtrll-, 0.5, POTNTS) 
CALL PLO* 00,5, 0.5, POINTS 7 

CALL PLOX 00.5, AH, 5), 0.5, POINTS) 
CALL PLOV 00.5, AftO), 0.5, POINTS) 
CALL PLOX (-0.5, AO,p), 0.5, POINTS) 




